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I.  INTRODUCTION 


Cavity  decoupled  underground  nuclear  explosions  in 
the  yield  range  from  1  to  10  kt  can  be  expected  to  generate 
seismic  signals  corresponding  to  m^  values  in  the  2.0  to  3.0 
range  and,  consequently,  such  events  pose  a  severe  test  for 
a  seismic  monitoring  network.  In  particular,  even  assuming 
the  existence  of  an  internal  network  capable  of  detecting 
such  events,  it  is  not  yet  clear  how  they  could  be  differen¬ 
tiated  from  the  background  of  many  small  earthquakes  and 
routine  chemical  explosions  of  comparable  size  which  can  be 
expected  to  be  detected  by  the  network.  Evernden  et  al. 

(1986)  and  others  have  suggested  that  it  may  be  possible  to 
resolve  this  discrimination  problem  using  high  frequency 
(>  20  Hz)  seismic  data.  However,  these  arguments  are  specu¬ 
lative  to  the  extent  that  little  is  currently  known  about  the 
high  frequency  source  characteristics  of  such  events,  parti¬ 
cularly  for  cavity  decoupled  explosions.  The  objectives  of 
the  research  summarized  in  this  report  have  been  to  conduct 
a  more  rigorous  theoretical  investigation  of  the  high  fre¬ 
quency  characteristics  of  decoupled  seismic  sources  and  to 
assess  the  implications  of  these  investigations  with  regard 
to  the  detection  and  discrimination  of  small  decoupled  explo¬ 
sions.  This  is  being  accomplished  through  a  combination  of 
detailed  deterministic  simulations  and  analyses  of  empirical 
data  recorded  from  selected  cavity  decoupled  explosions. 

The  organization  of  this  report  may  be  briefly  summa¬ 
rized  as  follows.  A  review  of  decoupling  phenomenology  is 
presented  in  Section  II,  with  specific  emphasis  on  those  as¬ 
pects  of  the  problem  affecting  high  frequency  coupling  effec¬ 
tiveness.  This  is  followed  in  Section  III  by  a  description 
of  the  results  of  a  series  of  detailed,  nonlinear  finite  dif¬ 
ference  calculations  which  have  been  carried  out  to  simulate 
the  seismic  source  functions  corresponding  to  coupled  and 
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decoupled  explosions  in  unsaturated  tuff  and  salt  emplacement 
media.  Results  of  preliminary  analyses  of  high  frequency 
seismic  data  recorded  from  the  STERLING  decoupled  explosion 
in  salt  and  the  MILL  YARD  decoupled  explosion  in  unsaturated 
tuff  are  presented  in  Section  IV.  This  is  followed  in  Section 
V  by  a  summary  and  a  statement  of  preliminary  conclusions  re¬ 
garding  high  frequency  seismic  monitoring  of  small  decoupled 
explosions. 


II.  REVIEW  OF  DECOUPLING  PHENOMENOLOGY 

The  concept  that  it  might  be  possible  to  significantly 
reduce  the  seismic  signal  radiated  by  an  underground  nuclear 
explosion  by  detonating  the  device  in  a  cavity  (i.e.,  cavity 
decoupling)  was  first  proposed  publicly  by  A.  L.  Latter  at 
the  1959  Nuclear  Test  Ban  Conference  in  Geneva  (Latter,  1959). 
However,  despite  the  fact  that  nearly  30  years  have  now  elapsed  . 
since  its  introduction,  a  number  of  major  issues  of  importance 
with  respect  to  seismic  monitoring  still  remain  unresolved. 

This  is  particularly  true  with  regard  to  the  potential  utility 
of  high  frequency  seismic  data  of  the  type  which  might  even¬ 
tually  be  collected  by  regional,  in-country  monitoring  sta¬ 
tions.  Before  proceeding  on  to  a  detailed  theoretical  consid¬ 
eration  of  these  issues,  it  is  appropriate  to  first  review 
some  of  the  general  principles  governing  the  seismic  source 
characteristics  of  coupled  and  decoupled  underground  nuclear 
explosions. 

Surrounding  every  fully  coupled  underground  explosion 
is  a  region  in  which  the  material  response  is  nonlinear. 

This  corresponds  to  the  regime  where  the  shock  pressure  is 
high  enough  to  cause  vaporization,  crushing  and  cracking  of 
the  medium.  As  the  range  from  the  detonation  point  increases, 
however,  the  shock  pressure  decays  to  a  level  (Pel)  at  which 
the  response  is  linear  and  it  is  the  forcing  function  acting 
at  this  radius  which  defines  the  characteristics  of  the  radi¬ 
ated  seismic  waves  which  are  used  for  detection  and  discrim¬ 
ination  purposes.  This  radius  is  commonly  called  the  elastic 

radius  and  denoted  as  r  , .  Consider  the  case  in  which  the 

el 

pressure  function  acting  at  rgl  is  a  simple  step  in  pressure, 

Pel  H  (t) .  Then,  it  can  be  shown  (Mueller  and  Murphy,  1971) 
that  the  amplitude  of  the  seismic  source  function,  S(w),  in 
the  low  frequency  limit  is  given  by 


S(fo) 


—  P  ,  r  , ' 
el  el 


Now,  if  Pe^  is  proportional  to  the  overburden  pressure,  pgh, 
then  it  follows  from  (1)  that,  at  a  fixed  depth  h,  the  low 
frequency  coupling  efficiency  will  depend  only  on  r  But, 
to  a  first  approximation,  the  peak  shock  pressure  decreases 
with  distance  r  as  r~n  where  n  is  the  medium-dependent  atten¬ 
uation  rate.  Thus,  rg^  decreases  as  n  increases,  and  the  rea¬ 
son  cavity  decoupling  works  is  that  for  strong  shocks  in  air 
n  =  3,  while  in  rock  n  «  2  (Brode,  1968;  Rodean,  1971).  That 
is,  am  explosion  is  said  to  be  fully  decoupled  if  it  is  det¬ 
onated  in  a  cavity  which  is  large  enough  that  the  surrounding 
medium  undergoes  no  nonlinear  deformation.  It  folxows  from 
the  above  discussion  that  the  radius  of  the  cavity  required 
to  fully  decouple  an  explosion  is  smaller  than  the  elastic 
radius  associated  with  a  fully  coupled  explosion  of  the  same 
yield  and  consequently,  by  equation  (1) ,  its  low  frequency 
coupling  efficiency  is  lower. 

The  criterion  for  full  decoupling  is  usually  ex¬ 
pressed  in  terms  of  a  requirement  that  the  late-time,  equilib¬ 
rium  pressure  in  the  cavity  be  less  than  or  equal  to  some  con¬ 
stant,  k,  times  the  overburden  pressure.  That  is  (Latter, 
et  al. ,  1961) 


(y-l)W 


i  it  r3 
3  c 


<  k  pgh 


where  y  is  the  gas  constant  (1.2-  for  air),  W  is  the  yield, 
rc  is  the  cavity  radius  and  k  is  generally  taken  to  lie  in 
the  range  from  0.5  (Latter  criterion)  to  1.0  (Patterson  cri¬ 
terion).  Thus,  the  Latter  criterion  incorporates  a  safety 
factor  of  two  which  was  designed  to  accommodate  the  large 
initial  pressure  spike  which  acts  on  the  cavity  wall  before 
the  cavity  pressure  reaches  the  equilibrium  value  given  by 
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equation  (2).  We  will  return  to  consider  the  effects  of  this 
spike  in  more  detail  later  in  this  report.  In  any  case,  it 
follows  from  (2)  that,  using  this  criterion,  the  minimum  cav¬ 
ity  radius  required  for  full  decoupling  will  be  given  by 


2.74  x  10' 


(k  p) 


for  W  in  kt,  h  in  m  and  p  in  gm/cm3 .  Thus,  the  cavity  radius 
required  for  full  decoupling  is  proportional  to  the  cube-root 
of  the  yield  and  inversely  proportional  to  the  cuke-root  of 
the  depth  of  burial.  This  is  illustrated  in  Figure  1  which 


shows  r  as  a  function  of  W  and  h  for  the  case  k  =  1,  p  = 
c  ' 

2.0  gm/cm?  The  depth  range  considered  extends  from  what  is 


denoted  "minimum  containment  depth"  to  a  fixed  (i.e.,  yield 
independent)  upper  limit  of  2000  m  which  corresponds  to  a 
constraint  (due  to  plastic  flow)  on  the  depth  at  which  stable 
cavities  can  be  constructed  in  salt.  The  minimum  containment 


depth  refers  to  the  depth  required  for  containment  of  a  coupled 
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explosion  of  yield  W  (i.e.,  approximately  122  W  m) .  Although 
this  criterion  does  not  apply  specifically  to  decoupled  explo¬ 
sions,  it  is  probably  representative  of  the  type  of  conserv¬ 
atism  an  evader  would  have  to  adopt.  The  dashed  horizontal 
lines  on  this  figure  denote  the  maximum  feasible  cavity  radii 
as  estimated  by  Evernden  (1976)  and  Rodean  (1971).  Evernden's 
upper  limit  is  based  primarily  on  theoretical  arguments  and 
appears  to  be  overly  pessimistic  given  that  the  Russians  have 
reported  the  existence  of  a  free-standing,  explosively- 
generated  cavity  in  salt  which  is  bigger  than  his  upper  bound 
(35  m  vs.  30  m) .  Rodean's  upper  limit  estimate,  on  the  other 
hand,  is  empirically  based  and  would  seem  to  be  more  nearly 
appropriate  in  that  it  corresponds  to  the  largest  cavities 
which  have  evt;r  been  constructed  by  the  natural  gas  industry 
using  solution  mining  techniques  in  salt.  However,  it  should 
be  noted  that  these  very  large  solution  cavities  are  generally 
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Figure  1.  Cavity  radii  required  for  full  decoupling  as  a 
function  of  source  depth  and  yield;  Patterson 
criteria,  p  =  2.0  gm/cm3 . 


neither  spherical  or  free-standing  (i.e.,  they  usually  are 
filled  with  pressurized  gas).  Thus,  the  feasibility  of  con¬ 
structing  decoupling  cavities  with  sizes  near  Rodean's  upper 
limit  has  yet  to  be  fully  demonstrated.  Moreover,  for  media 
other  than  salt,  there  is  no  evidence  that  free-standing  cav¬ 
ities  with  radii  approaching  even  Evernden's  limit  have  ever 
been  successfully  constructed.  Thus,  there  are  significant 
practical  limitations  on  the  size  of  an  explosion  which  can 
be  fully  decoupled  in  an  underground  cavity. 

The  simplified  coupling  analysis  represented  by 
equation  (1)  is  a  low  frequency  approximation  and  is  not  gen¬ 
erally  applicable  to  the  high  frequency  seismic  signals  which 
might  be  detected  by  a  regional  network.  The  reason  for  this 
is  illustrated  in  Figure  2  which  compares  estimates  of  the 
seismic  source  functions  for  coupled  and  fully  decoupled  10 
kt  explosions  at  a  depth  of  828  m  in  salt  (Murphy,  1979). 

It  can  be  seen  that,  using  these  analytic  approximations, 

both  source  functions  are  approximately  flat  at  low  frequen- 
.  -2 

cies  and  decrease  as  cu  above  a  characteristic  corner  fre¬ 
quency.  Moreover,  since  this  characteristic  corner  frequency 
is  proportional  to  rfil  1,  it  is  larger  for  the  decoupled  ex¬ 
plosion  than  it  is  for  the  coupled  explosion.  Thus,  taking 
the  ratio  of  these  two  source  functions  gives  rise  to  a  fre¬ 
quency  dependent  decoupling  factor  such  as  that  shown  in 
Figure  3.  It  can  be  seen  that  for  this  source  approximation 
the  decoupling  factor  is  approximately  constant  up  to  the  cor¬ 
ner  frequency  of  the  coupled  source  function,  above  which  it 

“2  * 
decreases  as  w  up  to  the  corner  frequency  associated  with 

the  decoupled  source  function.  Beyond  this  frequency,  the 

decoupling  remains  at  a  constant  (lower)  value,  or  at  least 

begins  to  decrease  less  rapidly,  depending  on  the  details  of 

the  high  frequency  specification  of  the  decoupled  source 

function.  This  prediction  that  decoupling  efficiency  will 

decrease  with  increasing  frequency  constitutes  one  of  the 
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reasons  why  high  frequency  monitoring  of  such  explosions  ap¬ 
pears  to  be  attractive.  On  the  other  hand,  as  has  been 
pointed  out  by  Murphy  (1979)  and  Evernden  et  ad.  (1986),  it 
may  be  feasible  for  an  evader  to  at  least  partially  overcome 
this  limitation  and  increase  the  decoupling  efficiency  at 
high  frequencies.  The  reasoning  that  leads  to  this  conclu¬ 
sion  may  be  summarized  as  follows.  Assume  that  the  decoupled 
seismic  source  function  can  be  approximated  as  a  simple  step 
in  pressure,  P  H(t),  acting  on  the  wall  of  a  cavity  which  is 
large  enough  to  fully  decouple  the  explosion.  Then,  by  anal¬ 
ogy  with  equation  (1) ,  it  can  be  shown  (Mueller  and  Murphy, 
1971)  that  the  amplitude  levels  of  the  seismic  source  func¬ 
tion  in  the  low  and  high  frequency  limits  scale  with  the 
pressure  level  and  cavity  radius  according  to  the  relations 


lim  |  S  (co)  |  ~  P  r' 


lim 

CO-voo 


(4) 


S  (co )  ~  P  r 


It  follows  that  if  r  is  increased  by  a  factor  k,  then  r 
increases  by  a  factor  k* .  However,  by  equation  (2),  the 
equilibrium  pressure  in  the  cavity  P  ~  r-3,  so  that  P  de- 
creases  by  ks  and  the  product  P  r^  is  independent  of  rc. 

Thus,  as  was  originally  concluded  by  Latter  et:  al.  (1961), 
the  low  frequency  decoupling  effectiveness  is  independent  of 
cavity  radius  as  long  as  rc  is  greater  than  the  radius  re¬ 
quired  for  full  decoupling,  rd.  However,  by  the  same  line 
of  reasoning,  it  can  be  seen  from  the  second  of  equations 
(4)  that  if  rc=  k  r^  (k  >  1) ,  then  the  high  frequency  de¬ 
coupled  seismic  source  decreases  by  a  factor  of  k2  and,  con¬ 
sequently,  the  high  frequency  decoupling  efficiency  increases 
by  this  amount.  This  predicted  effect  is  illustrated  in 
Figure  4  which  shows  the  frequency  dependent  decoupling 


expected  as  a  function  of  cavity  radius  for  a  10  kt  explosion 
in  salt.  Thus,  the  decoupling  efficiency  at  high  frequencies 
predicted  by  this  simple  source  model  can  be  improved  some¬ 
what  at  the  cost  of  constructing  a  larger  cavity. 

In  the  preceding  discussion  we  have  briefly  summa¬ 
rized  the  conventional  wisdom  concerning  the  seismic  source 
characteristics  of  cavity  decoupled  explosions  as  represented 
by  the  analytic  solution  to  the  idealized  problem  of  the 
linear,  elastic  response  of  an  infinite  homogeneous  medium  to 
a  step  in  pressure  in  a  spherical  cavity.  The  applicability 
of  this  model  at  low  frequencies  has  been  verified  to  some 
extent  by  the  limited  experimental  data  which  are  available 
for  analysis.  However,  it  has  long  been  recognized  that  this 
simple,  low  frequency  approximation  may  not  be  applicable  to 
the  high  frequencies  of  potential  interest  in  regional  seis¬ 
mic  monitoring.  For  example.  Figure  5  shows  a  comparison  of 
the  cavity  pressure  history  computed  by  Patterson  (1964)  for 
a  5  kt  explosion  in  a  51.2  m  radius  cavity  using  a  nonlinear 
finite  difference  code  with  the  simple  step  pressure  approxi¬ 
mation.  It  can  be  seen  that  the  more  precise  finite  differ¬ 
ence  solution  predicts  a  large  (nearly  1.6  kilobar)  initial 
pressure  spike  on  the  cavity  wall,  followed  by  a  rapid  decay 
and  subsequent  oscillation  about  the  steady-state  value  char¬ 
acteristic  of  the  step  function  approximation.  It  is  not 
clear  at  the  present  time  how  this  complex  pressure  history 
might  affect  the  corresponding  decoupled  seismic  source  func¬ 
tions,  but  Figure  6  shows  a  comparison  of  the  seismic  source 
functions  computed  from  the  two  pressure  histories  of  Figure 
5  assuming  linear,  elastic  response  of  the  medium  in  both 
cases  (Murphy,  1979).  It  can  be  seen  that  these  two  source 
estimates  are  roughly  comparable  at  low  frequency,  but  diverge 
significantly  for  frequencies  above  about  20  Hz.  For  exam¬ 
ple,  the  peak  at  about  65  Hz  in  Patterson's  (1964)  finite 
difference  source  corresponds  to  the  inverse  of  the  shock 
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6.  Comparison  of  decoupled  seismic  source  approxima¬ 
tions  for  a  5  kt  explosion  in  a  51.2  m  radius 
cavity  in  SALMON  salt. 
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reverberation  time  in  the  cavity  which  is  not  accounted  for 
in  the  step  approximation.  A  more  fundamental  concern  is 
that  the  large  initial  pressure  spike  may  well  induce  non¬ 
linear  response  in  the  surrounding  medium  which  could  signi¬ 
ficantly  modify  the  solution  at  both  high  and  low  frequencies. 
This  complex  issue  will  be  addressed  in  more  detail  in  the 
following  section  using  nonlinear  finite  difference  simula¬ 
tions  of  a  variety  of  decoupled  seismic  sources. 
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III.  NONLINEAR  FINITE  DIFFERENCE  SIMULATIONS  OF  DECOUPLED 
EXPLOSIONS  IN  SALT  AND  TUFF  EMPLACEMENT  MEDIA 


It  was  noted  in  the  preceding  section  that  most 
previous  analyses  of  cavity  decoupling  have  employed  the 
highly  oversimplified  step  pressure  approximation  to  describe 
the  decoupled  seismic  source  function.  The  objective  of  the 
analyses  summarized  in  this  section  has  been  to  utilize  the 
results  of  detailed,  nonlinear  finite  difference  simulations 
of  cavity  decoupled  explosions  to  better  define  the  broadband 
characteristics  of  the  seismic  signals  actually  generated  by 
such  sources.  A  series  of  numerical  simulations  of  decoupled 
explosions  in  unsaturated  tuff  and  salt  emplacement  media 
have  been  conducted  in  an  attempt  to  accomplish  this  objec¬ 
tive.  The  experimental  reference  point  for  the  tuff  analysis 
has  been  the  MILL  YARD  test  conducted  in  a  hemispherical  cav¬ 
ity  at  NTS.  For  salt,  the  experimental  reference  point  has 
been  the  STERLING  test  which  was  conducted  in  Mississippi  in 
the  cavity  created  by  the  SALMON  explosion.  In  both  cases, 
the  zone  size  in  the  finite  difference  calculations  was  se¬ 
lected  to  be  small  enough  that  frequency  components  well  above 
the  corner  frequencies  of  the  sources  were  accurately  repre¬ 
sented  in  the  solution.  That  is,  the  simulated  seismic  source 
functions  are  relevant  to  the  assessment  of  high  frequency 
seismic  monitoring  of  cavity  decoupled  tests. 

The  MILL  YARD  experiment  was  conducted  at  NTS  in  an 
11  m  radius  hemispherical  cavity  excavated  in  unsaturated 
tuff  at  a  depth  of  375  m.  The  overburden  pressure  at  this 
depth  is  about  65  bars  and,  consequently,  an  11  m  radius 
spherical  cavity  at  this  depth  would  be  expected  to  fully  de¬ 
couple  an  explosion  in  the  20  (Latter  criterion)  to  40 
(Patterson  criterion)  ton  yield  range.  Consequently,  non¬ 
linear  numerical  simulations  of  the  spherical  cavity  problem 
have  been  conducted  for  both  the  20  ton  and  40  ton  cases 
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using  finite  difference  zoning  fine  enough  to  accurately  rep¬ 
resent  frequency  components  as  high  as  500  Hz.  The  constitu¬ 
tive  model  for  tuff  used  in  these  nonlinear  simulations  was 
based  on  the  results  of  laboratory  measurements  on  core  sam¬ 
ples  taken  from  the  MILL  YARD  test  medium.  Shear  failure  was 
modeled  using  a  nonassociated  flow  rule  with  a  parabolic 
failure  surface  having  a  stress  difference  at  zero  mean  stress 
of  0.06  kilobars,  an  unconfined  compressive  strength  of  0.10 
kilobars  and  a  maximum  stress  difference  of  0.25  kilobars  at 
mean  stress  levels  above  0.30  kilobars.  Based  on  the  results 
of  in-situ  sonic  logging,  the  P  and  S  wave  velocities  were 
taken  to  be  2558  and  1372  m/sec,  respectively.  A  detailed 
air  equation  of  state  was  used  to  model  the  shock  propagation 
inside  of  the  11  m  radius  cavity. 

The  cavity  pressure  as  a  function  of  time  computed 
at  a  grid  point  near  the  cavity  wall  for  the  20  ton  nonlinear 
simulation  is  shown  in  Figure  7.  It  can  be  seen  that  this 
computed  pressure  profile  is  similar  to  that  shown  previously 
in  Section  II  (cf.  Figure  5)  in  that  it  is  characterized  by  a 
sharp  initial  pressure  spike  associated  with  the  reflection 
of  the  shock  from  the  cavity  wall,  followed  by  reverberation 
oscillations  which  gradually  damp  out  as  the  pressure  converges 
to  the  late-time  equilibrium  value  given  by  equation  (2).  The 
Fourier  transform  of  this  pressure  loading  is  shown  in  Figure 
8  where  it  is  compared  with  that  corresponding  to  the  simple 
step  pressure  approximation.  As.  was  noted  previously,  the 
step  function  approximation  significantly  underestimates  the 
high  frequency  content  of  the  computed  pressure  acting  on  the 
cavity  wall,  particularly  near  the  frequencies  corresponding 
to  the  reverberation  of  the  shock  in  the  cavity. 

The  seismic  source  function  computed  for  the  20  ton 
nonlinear  simulation  is  shown  in  Figure  9  where  it  is  com¬ 
pared  with  the  linear  elastic  solution  corresponding  to  the 
same  cavity  pressure  loading.  It  can  be  seen  that  these  two 
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Figure  8.  Comparison  of  the  spectral  composition  of  the  com¬ 
plex  pressure  profile  of  Figure  7  with  that  asso¬ 
ciated  with  the  simple  step  pressure  approximation 
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Figure  9.  Comparison  of  linear  and  nonlinear  seismic  source 
functions  computed  for  a  20  ton  decoupled  explo¬ 
sion  in  an  11  m  radius  spherical  cavity  in  un¬ 
saturated  tuff. 


source  functions  are  generally  quite  similar,  particularly 
at  high  frequencies,  with  both  showing  comer  frequencies  of 
about  35  Hz  and  strong  secondary  reverberation  peaks  at  around 
200  Hz.  Somewhat  surprisingly,  however,  this  comparison  in¬ 
dicates  increased  decoupling  effectiveness  at  low  frequency 
in  the  nonlinear  solution  over  that  predicted  for  the  linear, 
elastic  solution,  in  apparent  contradiction  of  the  decoupling 
concept.  Detailed  examination  of  the  nonlinear  solution  in¬ 
dicates  that  this  phenomenon  can  be  explained  by  the  fact  that 
in  this  unsaturated  tuff  medium,  crush-up  of  the  air  filled 
porosity  (1.8  percent)  results  in  a  decrease  in  the  computed 
permanent  displacement  at  the  elastic  radius  and  a  correspond¬ 
ing  reduction  in  the  low  frequency  coupling  effectiveness. 

In  an  attempt  to  gain  more  insight  into  the  specific  cause  of 
this  reduction,  supplemental  20  ton  simulations  were  performed 
in  which  the  seismic  source  functions  corresponding  to  the 
classical  step  pressure  approximation  were  estimated  using 
both  the  nonlinear  and  linear  elastic  models  of  the  tuff  re¬ 
sponse.  In  this  case,  the  two  computed  source  functions  were 
found  to  be  essentially  identical  and  it  was  therefore  con¬ 
cluded  that  the  nonlinear  reduction  in  low  frequency  seismic 
coupling  noted  in  Figure  9  must  be  associated  with  the  re¬ 
sponse  to  the  initial  pressure  spike  acting  on  the  cavity 
wall.  It  follows  that  the  safety  factor  in  the  Latter  de¬ 
coupling  criterion,  which  was  designed  to  accommodate  this 
pressure  spike,  is  not  sufficient  to  insure  linear,  elastic 
response,  at  least  in  unsaturated,  porous  tuff  media. 

The  seismic  source  function  computed  for  a  40  ton 
explosion  in  the  same  11  m  radius  spherical  cavity  in  tuff 
using  the  nonlinear  finite  difference  algorithm  is  shown  in 
Figure  10  where  it  is  compared  with  the  corresponding  20  ton 
solution  which  has  been  multiplied  by  a  factor  of  two  for 
approximate  yield  normalization.  This  comparison  indicates 
that  the  relative  coupling  effectiveness  at  40  tons  is  only 
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Comparison  of  40  ton  and  yield-normalized  20 
ton  nonlinear  seismic  source  functions  computed 
for  decoupled  explosions  in  an  11  m  radius 
spherical  cavity  in  unsaturated  tuff. 


slightly  increased  at  low  frequency.  In  fact,  the  only  signi¬ 
ficant  difference  between  the  two  solutions  is  in  the  fre¬ 
quency  of  the  reverberation  peak  which  shifts  to  higher  fre¬ 
quency  in  the  40  ton  case  due  to  the  larger  shock  velocities. 
These  results  suggest  that  there  is  no  reason  to  prefer  the 
more  conservative  Latter  decoupling  criterion  over  that  pro¬ 
posed  by  Patterson,  at  least  for  this  source  medium.  The  fre¬ 
quency  dependent  decoupling  factors  computed  from  these  20  and 
40  ton  simulations  are  shown  in  Figure  11  where  it  can  be  seen 
that,  as  expected,  they  are  nearly  constant  at  low  frequency 
and  decrease  with  increasing  frequency  above  the  corner  fre¬ 
quencies  associated  with  the  corresponding  coupled  source 
functions.  The  principal  differences  with  respect  to  the  sim¬ 
ple  step  function  approximation  are  the  pronounced  minima  in 
the  decoupling  effectiveness  around  200  Hz  associated  with  the 
cavity  reverberation  peaks.  The  computed  decoupling  levels 
decrease  from  about  a  factor  of  40  at  low  frequency  to  values 
less  than  one  above  about  200  Hz.  These  values  are  low  rela¬ 
tive  to  those  obtained  from  comparable  simulations  in  salt 
and  reflect  the  rapid  shock  attenuation  which  is  predicted  to 
occur  in  this  weak  tuff  medium  under  the  high  pressures  in¬ 
duced  by  coupled  explosions.  However,  with  respect  to  detect¬ 
ability,  it  is  appropriate  to  note  that  the  m^  values  for 
events  in  dry  tuff  have  been  observed  to  be  as  much  as  a  full 
order  of  magnitude  lower  than  those  associated  with  explosions 
of  comparable  yield  in  hardrock  or  saturated  media.  Thus, 
very  low  signal  levels  would  be  expected  from  small  decoupled 
explosions  in  this  medium. 

The  STERLING  experiment  was  conducted  in  a  salt  dome 
in  Mississippi  in  the  cavity  produced  by  the  SALMON  explosion. 
The  cavity  was  approximately  spherical,  with  a  radius  of  about 
17  m  and  was  located  at  a  depth  of  828  m.  Since  the  overbur¬ 
den  pressure  at  this  depth  is  about  180  bars,  such  a  17  m 
radius  spherical  cavity  would  be  expected  to  fully  decouple 
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an  explosion  in  the  210  (Latter  criterion)  to  420  (Patterson 
criterion)  ton  yield  range.  STERLING  had  a  yield  of  380  tons, 
and  thus  was  close  to  the  Patterson  limit.  Consequently,  a 
nonlinear  simulation  of  a  380  ton  nuclear  explosion  in  a  17  m 
radius  spherical  cavity  in  salt  was  conducted  in  a  preliminary 
attempt  to  investigate  the  high  frequency  seismic  source  char¬ 
acteristics  of  decoupled  explosions  in  this  medium.  As  in  the 
tuff  simulations,  these  calculations  were  carried  out  using  a 
finite  difference  zoning  which  was  fine  enough  to  accurately 
represent  frequency  components  as  high  as  500  Hz.  The  consti¬ 
tutive  model  used  in  these  calculations  was  primarily  based 
on  the  results  of  laboratory  material  property  tests  on  poly¬ 
crystalline  salt  (Heard  et  al. ,  1975).  Shear  failure  was 
modeled  using  a  nonassociated  flow  rule  with  a  parabolic  fail¬ 
ure  surface  having  a  stress  difference  at  zero  mean  stress  of 
0.14  kilobars,  an  unconfined  strength  of  0.31  kilobars  and  a 
maximum  stress  difference  of  0.68  kilobars  at  mean  stress 
levels  above  0.60  kilobars.  P  and  S  wave  velocities  of  4551 
and  2440  m/sec,  respectively,  were  selected  on  the  basis  of 
analyses  of  SALMON  free-field  observations.  In  addition,  a 
detailed  equation  of  state  was  used  to  model  the  shock  propa¬ 
gation  inside  of  the  17  m  radius  cavity. 

The  cavity  pressure  as  a  function  of  time  computed 
at  a  grid  point  near  the  cavity  wall  for  the  380  ton  explo¬ 
sion  simulation  is  shown  in  Figure  12  where  it  can  be  seen 
that  it  is  qualitatively  similar  to  that  estimated  for  the 
tuff  simulation  (cf.  Figure  7).  However,  in  this  case  the 
initial  pressure  spike  has  an  amplitude  of  about  4  kilobars 
in  contrast  to  the  400  bar  spike  amplitude  of  Figure  7.  This 
reflects  the  fact  that  STERLING  was  conducted  at  much  greater 
depth  than  MILL  YARD  and,  therefore,  could  accommodate  signi¬ 
ficantly  higher  pressures  under  the  standard  decoupling  cri¬ 
teria.  That  is,  according  to  these  criteria,  a  significantly 
larger  (i.e.,  23  m  radius)  cavity  would  be  required  to  decouple 


380  tons  at  MILL  YARD  depth  and,  in  that  case,  much  lower 
peak  pressures  would  be  expected. 

The  seismic  source  function  computed  from  the  non¬ 
linear  simulation  is  shown  in  Figure  13  where  it  is  compared 
with  the  linear  elastic  solution  corresponding  to  the  same 
cavity  pressure  loading.  It  can  be  seen  that,  as  with  the 
tuff  simulations,  these  two  source  functions  are  quite  similar 
at  high  frequencies  with  both  showing  corner  frequencies  of 
about  40  Hz  and  strong  secondary  reverberation  peaks  at  around 
300  Hz.  Unlike  the  tuff  simulations,  however,  the  low  fre¬ 
quency  differences  in  this  case  are  in  accord  with  what  would 
be  expected  from  simple  decoupling  theory.  That  is,  the  non¬ 
linear  solution  predicts  slightly  (**  20  percent)  greater 
coupling  than  that  computed  using  the  linear  approximation. 
Once  again,  the  nonlinear  response  is  associated  with  the 
initial  pressure  spike  on  the  cavity  wall,  as  evidenced  by 
the  fact  that  the  computed  nonlinear  and  linear  responses  to 
the  classical  160  bar  step  pressure  loading  are  essentially 
identical.  Analysis  of  the  nonlinear  simulation  results  in¬ 
dicates  that  the  predominant  nonlinear  mechanism  in  this  case 
was  yielding  of  the  salt  out  to  a  radius  of  about  19  m,  as 
opposed  to  the  pore  crush-up  mechanism  which  predominated  in 
the  unsaturated  tuff  simulations.  Furthermore,  comparing 
Figures  13  and  9  it  can  be  seen  that  the  effects  of  the  non¬ 
linearity  on  the  seismic  coupling  are  much  less  pronounced  in 
these  salt  simulations  and,  in  fact,  simulations  conducted 
with  yields  closer  to  the  210  ton  limit  suggested  by  the 
Latter  criterion  indicate  little  or  no  nonlinear  response  of 
the  cavity  wall.  Thus,  for  this  salt  model,  there  is  a  dif¬ 
ference  between  the  decoupling  effectiveness  obtained  using 
the  Latter  versus  the  Patterson  decoupling  criterion.  How¬ 
ever,  the  differences  are  confined  to  low  frequency  and  do 
not  appear  to  be  very  significant  in  this  case. 


functions  computed  for  a  380  ton  decoupled  explo¬ 
sion  in  a  17  m  radius  spherical  cavity  in  salt. 


In  summary,  preliminary  cavity  decoupling  simula¬ 
tions  have  been  conducted  for  explosions  in  unsaturated  tuff 
and  salt  emplacement  media.  These  simulations  have  confirmed 
the  fact  that  the  simple  step  pressure  approximation  is  not 
valid  at  high  frequencies  and  that  the  initial  pressure  spike 
on  the  cavity  wall  can  induce  significant  nonlinear  response. 
However,  in  all  the  simulations  conducted  to  date  the  use  of 
the  more  realistic  cavity  pressure  loading  has  resulted  in 
significantly  enhanced  seismic  coupling  at  high  frequency  over 
that  predicted  by  the  simple  step  pressure  approximation.  It 
now  remains  to  be  determined  whether  these  predicted  high  fre¬ 
quency  effects  can  be  detected  in  the  observed  seismic  data 
and  used  for  identification  purposes. 


IV.  HIGH  FREQUENCY  CHARACTERISTICS  OF  THE  STERLING  AND 
MILL  YARD  SEISMIC  SOURCE  FUNCTIONS 

The  results  of  the  finite  difference  simulations 
which  were  described  in  the  preceding  section  indicated  that 
there  may  be  some  characteristics  of  the  high  frequency  seis¬ 
mic  source  functions  associated  with  cavity  decoupled  explo¬ 
sions  which  could  be  exploited  for  identification  purposes. 

In  this  section,  an  attempt  will  be  made  to  assess  whether 
evidence  of  these  potentially  diagnostic  characteristics  can 
be  identified  and  confidently  extracted  from  the  available 
experimental  data.  More  specifically,  near-field  seismic 
data  recorded  from  the  STERLING  and  MILL  YARD  nuclear  de¬ 
coupling  tests  will  be  examined  for  evidence  of  these  theo¬ 
retically  predicted  high  frequency  source  characteristics. 

4.1  STERLING 

STERLING  was  a  0.38  kt  nuclear  test  conducted  in 
December,  1966  in  the  cavity  created  by  the  SALMON  explosion 
which  had  previously  been  detonated  at  a  depth  of  828  m  in 
the  Tatum  salt  dome  near  Hattiesburg,  Mississippi.  As  is 
indicated  in  Figure  14,  the  cavity  was  approximately  spheri¬ 
cal  with  a  radius  of  about  17  m.  With  reference  to  equation 
(2),  0.38  kt  in  a  17  m  radius  spherical  cavity  would  be  ex¬ 
pected  to  result  in  a  late-time  equilibrium  pressure  level  of 
about  160  bars,  which  is  close  to  the  value  of  the  overburden 
pressure  of  180  bars  at  STERLING  depth.  Thus,  STERLING  was 
approximately  consistent  with  the  Patterson  decoupling  crite¬ 
rion. 

Most  previous  analyses  of  STERLING  free-field  data 
have  focused  on  lower  frequency  recordings  made  at  ranges  of 
greater  them  160  m.  However,  Sisemore  et  al.  (1969)  did  re¬ 
port  broadband  results  obtained  from  an  accelerometer  located 
near  shot  depth  at  a  range  of  only  52  m.  The  radial  component 
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velocity  time  history  inferred  from  this  recording  is  shown 
in  Figure  15  where  it  is  compared  with  the  results  of  two 
different  theoretical  simulations.  The  simulated  velocity 
profile  shown  in  the  middle  of  this  figure  corresponds  to  the 
computed  linear  elastic  response  of  the  salt  emplacement 
medium  to  the  complex  cavity  pressure  loading  determined  from 
a  spherically  symmetric  finite  difference  calculation,  while 
that  shown  at  the  bottom  of  the  figure  corresponds  to  the 
linear  elastic  response  to  the  classical  step  pressure  ap¬ 
proximation.  It  can  be  seen  from  this  comparison  that  even 
ignoring  anelastic  attenuation  effects,  the  solution  corre¬ 
sponding  to  the  160  bar  step  pressure  underestimates  the  ob¬ 
served  peak  velocity  by  about  a  factor  of  three  and  does  not 
match  the  oscillatory  character  of  the  observed  velocity  time 
history.  The  solution  corresponding  to  the  finite  difference 
pressure  simulation,  on  the  other  hand,  provides  a  much  bet¬ 
ter  qualitative  description  of  the  observed  data,  although  it 
does  somewhat  overpredict  the  peak  amplitude  and  dominant 
frequency  of  the  initial  pulse.  This  is  not  surprising,  since 
the  combined  effects  of  nonlinear  and  linear  anelastic  atten¬ 
uation  occurring  between  the  cavity  wall  and  the  observation 
point  at  52  m  have  been  ignored  in  this  simulation,  and  in¬ 
clusion  of  any  such  effects  would  be  expected  to  both  broaden 
the  initial  pulse  and  reduce  its  amplitude.  In  any  case, 
this  observation  provides  convincing  evidence  that  the  high 
frequency  characteristics  of  the  actual  cavity  pressure  load¬ 
ing  were  effectively  coupled  into  the  seismic  regime  in  this 
case. 

4 . 2  MILL  YARD 

The  MILL  YARD  nuclear  test  was  conducted  at  NTS  on 
October  9,  1985  in  a  hemispherical  cavity  excavated  in  tuff. 
The  cavity  configuration  is  illustrated  in  Figure  16  where 
it  can  be  seen  that  the  cavity  had  a  radius  of  11  m  and 


overlay  a  floor  consisting  of  a  5  in  thick  layer  of  unconsoli¬ 
dated  soil  (Garbin,  1986).  The  cavity  was  constructed  at  a 
depth  of  375  m  in  an  unsaturated  tuff  medium  characterized 
by  a  compressional  wave  velocity  of  about  2.5  km/sec.  Figure 
17  shows  a  vertical  section  through  the  shot  point  which  in¬ 
dicates  that  the  propagation  path  between  the  cavity  and  the 
ground  surface  directly  above  the  shot  (i.e.,  ground  zero) 
consists  of  a  variety  of  tuff  layers  characterized  by  differ¬ 
ent  physical  properties.  Note  that  although  there  is  some 
evidence  of  lateral  heterogeneity,  this  section  does  not  ap¬ 
pear  to  be  particularly  complex,  at  least  at  the  scale  of 
this  figure.  We  will  return  to  this  point  later  in  our  analy¬ 
sis  of  the  characteristics  of  the  seismic  motions  recorded 
at  ground  zero  from  MILL  YARD. 

Garbin  (1986)  analyzed  both  free-field  and  surface 
data  recorded  from  this  event  and  noted  that  the  spectral 
compositions  of  the  ground  motions  recorded  in  these  two 
regimes  are  quite  different.  That  is,  the  free-field  data 
recorded  at  distances  of  about  6  and  13  m  from  the  cavity 
wall  are  dominated  by  very  high  frequency  (~  1000  Hz)  compo¬ 
nents,  while  those  recorded  at  ground  zero  at  a  range  of 
375  m  are  dominated  by  components  with  frequencies  of  less 
than  about  50  Hz.  Garbin  concluded  that  the  MILL  YARD  free- 
field  data  are  dominated  by  the  motions  induced  by  the  initial 
pressure  spike  on  the  cavity  wall.  That  is,  as  with  STERLING, 
there  is  evidence  that  the  potentially  diagnostic,  high  fre¬ 
quency  cavity  source  components  do  effectively  couple  into 
the  seismic  regime.  However,  the  lower  frequency  content  of 
the  surface  recordings  indicates  that  frequency  dependent 
attenuation  can  mask  these  source  characteristics,  even  over 
relatively  short  propagation  paths.  This  observation  moti¬ 
vated  a  more  detailed  analysis  of  the  spectral  composition  of 
the  MILL  YARD  surface  recordings. 


The  initial  P  wave  displacement  spectrum  derived 
from  the  vertical  component  ground  motion  recorded  at  ground 
zero  is  shown  in  Figure  18  where  it  is  compared  with  the  spec¬ 
trum  predicted  theoretically  for  a  step  in  pressure  in  an  11  m 
radius  spherical  cavity  in  a  medium  characterized  by  a  com- 
pressional  wave  velocity  of  2.5  km/sec.  It  can  be  seen  that 
while  there  is  evidence  of  some  potentially  interesting  high 
frequency  complexity  in  this  observed  spectrum,  the  most  no¬ 
table  features  of  this  comparison  are  the  facts  that  the  ob¬ 
served  corner  frequency  is  lower  than  expected  and  the  spec¬ 
tral  decay  above  the  corner  frequency  is  much  more  rapid  than 
that  predicted  by  the  simple  elastic  solution.  It  was  ori¬ 
ginally  thought  that  these  discrepancies  might  be  related  to 
the  hemispherical  geometry  of  the  MILL  YARD  cavity.  However, 
subsequent  investigation  has  provided  strong  evidence  that 
these  two  spectra  can,  in  fact,  be  reconciled  within  the  con¬ 
text  of  a  fairly  simple  propagation  path  model. 

The  first  step  in  the  analysis  was  to  estimate  and 
correct  for  the  effects  of  anelastic  attenuation.  This  was 
accomplished  using  a  model  in  which  it  was  assumed  that  the 
P  wave  displacement  source  function  should  decrease  approxi¬ 
mately  as  co~2  above  the  source  corner  frequency  and  that  the 
frequency  dependent  attenuation  effects  of  propagation  could 
be  accounted  for  by  a  conventional,  constant  Q  operator. 
Application  of  this  model  to  the  observed  P  wave  spectrum  of 
Figure  18  leads  to  an  estimated  average  Q  value  of  about  10 
for  the  path  between  the  cavity  surface  and  ground  zero. 

Figure  19  shows  the  comparison  with  results  from  applying 
this  attenuation  correction  to  the  observed  spectrum  of  Fig¬ 
ure  18.  it  can  be  seen  that  the  Q-normalized  observed  spec¬ 
trum  is  fairly  consistent  with  the  simple  theoretical  predic¬ 
tion,  although  there  remain  some  significant  spectral  modula¬ 
tion  effects  which  are  not  accounted  for  by  the  model.  Analy¬ 
sis  of  the  corresponding  time  domain  signals  provides  insight 


into  the  cause  of  this  spectral  complexity.  For  example. 
Figure  20  shows  comparisons  between  the  observed  ground  zero 
acceleration,  velocity  and  displacement  waveforms  with  the 
vertical  component  synthetics  predicted  using  the  simple 
spherical  cavity  source  model  and  an  anelastic  propagation 
model  characterized  by  a  Q  value  of  10.  It  can  be  seen  from 
this  figure  that  the  dominant  periods  of  the  first  half  cycles 
of  the  observed  acceleration  and  velocity  waveforms  are 
matched  quite  well  by  these  synthetics,  indicating  that  the 
inferred  combination  of  source  and  Q  models  is  approximately 
correct.  However,  the  displacement  comparison  is  less  satis¬ 
factory  and  examination  of  the  acceleration  and  velocity  com¬ 
parisons  suggests  that  this  may  be  due  to  interference  with 
a  secondary  phase  of  opposite  polarity  arriving  about  0.03 
seconds  after  the  first  motion.  The  onset  of  this  arrival  is 
indicated  by  the  vertical  arrows  above  these  traces  and  its 
phase  relationship  to  the  first  arrival  is  most  clearly  indi¬ 
cated  in  the  comparison  of  the  observed  and  synthetic  veloc¬ 
ity  waveforms.  These  observations  suggest  that  the  modulation 
evident  in  the  observed  displacement  spectrum  of  Figure  19 
may  be  due  to  interference  between  arrivals.  In  fact,  the 
observed  short  time  delay  and  phase  relationship  between  the 
arrivals  is  suggestive  of  a  near-surface  reverberation  effect. 
The  plausibility  of  this  interpretation  is  demonstrated  more 
clearly  in  Figure  21  where  the  result  of  dividing  the  attenu¬ 
ation  corrected  spectrum  of  Figure  19  by  the  source  spectrum 
corresponding  to  the  step  in  pressure  in  the  11  m  radius 
spherical  cavity  is  displayed  and  compared  with  theoretical 
layer  response  functions  computed  assuming  a  two-way  travel 
time  in  the  layer  consistent  with  the  observed  interference 
pattern  and  layer  Q  values  of  «  and  10  respectively.  Thus, 
for  example,  if  the  near-surface  compressional  wave  velocity 
is  1000  m/sec,  then  the  observed  delay  time  of  0.03  seconds 
would  imply  a  layer  thickness  of  about  15  m,  which  is  not  an 
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unreasonable  value.  It  can  be  seen  from  this  figure  that  the 
Q  =  10  theoretical  layer  response  agrees  quite  well  with  the 
normalized  observed  spectrum  below  about  80  Hz.  The  extent 
to  which  this  model  fits  the  time  domain  data  is  illustrated 
in  Figure  22  where  the  comparison  shown  previously  in  Figure 
20  has  been  modified  to  incorporate  the  effect  of  this  Q  =  10 
layer  response  into  the  synthetics.  It  can  be  seen  that  these 
initial  portions  of  the  acceleration,  velocity  and  displace¬ 
ment  synthetic  waveforms  are  now  in  excellent  agreement  with 
the  observations,  confirming  the  applicability  of  the  inferred 
site  response  model.  Finally,  Figure  23  shows  the  comparison 
between  the  attenuation-normalized  observed  spectrum  and  the 
theoretical  spectrum  obtained  by  multiplying  the  step  func¬ 
tion  solution  times  the  layer  response  function.  This  com¬ 
parison  seems  to  confirm  the  fact  that  the  MILL  YARD  P  wave 
source  function  can  be  adequately  approximated  by  a  simple 
step  in  pressure  in  an  11  m  radius  spherical  cavity,  at  least 
over  this  limited  frequency  range.  That  is,  as  was  suggested 
by  the  finite  difference  simulations  described  in  the  pre¬ 
vious  section  (cf .  Figure  10) ,  higher  frequency  data  will  be 
required  to  confidently  identify  the  effects  of  the  pressure 
spike  and  reverberations  in  the  MILL  YARD  cavity.  Unfortu¬ 
nately,  the  high  degree  of  anelastic  attenuation  occurring 
along  the  propagation  path  to  ground  zero  does  not  allow  us 
to  recover  this  high  frequency  information  in  this  case. 

It  follows  from  the  above  discussion  that  very  de¬ 
tailed  path  corrections  will  be  required  if  high  frequency 
seismic  data  are  to  be  used  to  identify  small  de'coupled  ex¬ 
plosions.  Another  example  of  the  complexity  of  such  propaga¬ 
tion  effects  is  provided  by  Figure  24  which  shows  the  verti¬ 
cal  and  two  horizontal  components  of  displacements,  velocity 
and  acceleration  recorded  at  ground  zero  from  MILL  YARD.  In 
each  case,  the  three  orthogonal  components  of  motion  are 
plotted  at  the  same  absolute  amplitude  scale  so  that  the 
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Figure  22. 


Comparison  of  observed  MILL  YARD  ground  zero 
acceleration,  velocity  and  displacement  waveforms 
with  the  vertical  component  synthetics  computed 
by  convolving  the  synthetics  of  Figure  20  with 
the  Q  =  10  layer  response  function  from  Figure 
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vector  character  of  the  motion  can  be  examined.  Now,  in  a 
plane-layered  medium  the  initial  P  wave  motion  at  ground  zero 
would  be  expected  to  be  negligible  on  the  horizontal  compo¬ 
nents,  and  it  can  be  seen  that  this  is  essentially  true  for 
the  displacement  time-histories.  However,  in  the  progression 
to  the  higher  frequency  velocity  and  acceleration  traces, 
this  approximation  begins  to  break  down  until,  for  the  accel¬ 
eration  motion,  the  initial  peak  actually  occurs  on  one  of 
the  horizontal  components.  This  is  rather  surprising  given 
the  relatively  uncomplicated  vertical  section  through  this 
propagation  path  shown  previously  in  Figure  17  and  provides 
further  illustration  of  some  of  the  complications  which  can 
obscure  the  interpretation  of  high  frequency  seismic  data  for 
identification  purposes. 
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V.  SUMMARY  AND  CONCLUSIONS 


5 . 1  SUMMARY 

Recently,  Evernden  et_  al_.  (19  86)  and  others  have 
suggested  that  the  detection  and  identification  of  cavity  de¬ 
coupled  nuclear  explosions  can  best  be  carried  out  by  seismic 
monitoring  at  frequencies  above  20  Hz.  However,  this  hypoth¬ 
esis  is  based,  at  least  in  part,  upon  the  predictions  of  a 
highly  oversimplified  analytic  seismic  source  model  for  de¬ 
coupled  explosions,  and  the  applicability  of  this  model  at 
high  frequencies  has  yet  to  be  demonstrated  either  theoreti¬ 
cally  or  experimentally.  The  investigations  summarized  in 
this  interim  report  have  centered  on  more  rigorous  investiga¬ 
tions  of  the  high  frequency  characteristics  of  such  decoupled 
seismic  sources.  More  specifically,  a  combination  of  detailed 
deterministic  simulations  and  analyses  of  empirical  data  re¬ 
corded  from  selected  cavity  decoupled  explosions  has  been 
employed  to  investigate  this  issue. 

A  brief  review  of  decoupling  phenomenology  was  pre¬ 
sented  in  Section  II  where  the  various  decoupling  criteria 
were  defined  and  used  to  illustrate  the  dependence  of  the  re¬ 
quired  cavity  size  on  explosion  yield  and  depth  of  burial. 

In  addition,  the  classical  step  function  pressure  approxima¬ 
tion  to  the  decoupled  seismic  source  function  was  compared 
with  the  seismic  source  corresponding  to  the  more  complex 
cavity  pressure  history  predicted  using  nonlinear  finite  dif¬ 
ference  calculations  in  order  to  illustrate  the  discrepancies 
between  these  two  source  approximations  at  high  frequencies. 

The  results  of  a  series  of  detailed,  nonlinear 
finite  difference  simulations  of  the  seismic  source  functions 
corresponding  to  cavity  decoupled  explosions  in  unsaturated 
tuff  and  salt  emplacement  media  were  described  in  Section  III 
where  they  were  compared  with  the  results  of  corresponding 
linear  simulations.  These  comparisons  were  then  used  to 
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identify  the  mechanisms  by  which  the  initial  transient  pres¬ 
sure  spike  on  the  cavity  wall  induces  nonlinear  responses  in 
these  two  media.  As  a  result  of  this  analysis,  it  was  demon¬ 
strated  that  the  use  of  a  more  realistic  cavity  pressure  load¬ 
ing  leads  to  significantly  enhanced  estimates  of  the  seismic 
coupling  at  high  frequency  relative  to  that  predicted  by  the 
simple  step  pressure  approximation. 

Preliminary  analyses  of  high  frequency  seismic  data 
recorded  from  the  STERLING  decoupled  explosion  in  salt  and 
the  MILL  YARD  decoupled  explosion  in  unsaturated  tuff  were 
presented  in  Section  IV  where  an  attempt  was  made  to  determine 
whether  the  theoretically  predicted  high  frequency  decoupling 
effects  described  in  Section  III  could  be  detected  in  the  ob¬ 
served  seismic  data.  Using  the  results  of  these  analyses,  it 
was  demonstrated  that,  while  the  high  frequency  source  compo¬ 
nents  associated  with  the  complex  cavity  pressure  loading  do 
effectively  couple  into  the  seismic  regime,  the  effects  of 
propagation  make  it  difficult  to  reliably  extract  these  fea¬ 
tures  from  the  measured  data. 

5.2  CONCLUSIONS 

The  results  of  the  analyses  described  above  support 
the  following  preliminary  conclusions  concerning  the  high  fre¬ 
quency  seismic  source  characteristics  of  cavity  decoupled 
underground  nuclear  explosions. 

(1)  The  results  of  a  series  of  preliminary  non¬ 
linear  finite  difference  simulations  of  cavity  decoupling  in 
tuff  and  salt  emplacement  media  have  confirmed  that  the  sim¬ 
ple  step  pressure  approximation  does  not  provide  an  adequate 
description  of  the  high  frequency  seismic  characteristics  of 
such  sources.  In  particular,  all  the  simulations  conducted 
to  date  indicate  enhanced  seismic  coupling  at  high  frequency 
relative  to  that  predicted  using  the  classical  step  pressure 
approximation . 
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(2)  The  nonlinear  simulations  also  indicate  that, 
while  the  Latter  decoupling  criterion  may  be  approximately 
correct  for  salt,  it  is  not  adequate  to  prevent  significant 
pressure-spike  induced  nonlinearities  in  unsaturated  tuff. 

For  both  tuff  and  salt,  however,  the  decoupling  effectiveness 
obtained  using  the  Patterson  criterion  is  not  significantly 
different  from  that  obtained  using  the  more  conservative 
Latter  criterion. 

(3)  The  nature  of  the  nonlinear  response  in  unsat¬ 
urated  tuff  is  such  that  the  decoupling  effectiveness  at  low 
frequency  is  actually  increased  over  that  predicted  by  the 
linear,  elastic  solution.  This  reflects  the  fact  that  the 
predominant  nonlinear  mechanism  in  this  medium  is  crush-up 

of  the  air  filled  porosity.  In  salt,  on  the  other  hand,  where 
the  predominant  nonlinear  mechanism  is  yielding,  the  nonlinear 
solution  predicts  decreased  decoupling  effectiveness  relative 
to  the  linear,  in  accord  with  traditional  decoupling  concepts. 

(4)  Examination  of  the  free-field  data  recorded 
from  both  the  STERLING  and  MILL  YARD  nuclear  decoupling  ex¬ 
periments  indicates  that  the  high  frequency  source  components 
associated  with  the  actual  complex  cavity  pressure  loadings 
do  effectively  couple  into  the  seismic  regime,  in  agreement 
with  the  predictions  of  the  nonlinear  finite  difference  simu¬ 
lations. 

(5)  However,  the  observed  complexity  of  the  propa¬ 
gation  effects  at  high  frequency,  as  evidenced  by  the  MILL 
YARD  ground  zero  data,  raises  questions  as  to  whether  such, 
source  characteristics  can  be  reliably  recovered  from  the 
recorded  seismic  data  and  used  for  identification  purposes. 
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